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The analytical and numerical study of the angular dependence of magnetoresistance in layered 
quasi-two-dimensional (Q2D) metals is performed. The harmonic expansion analytical formulas for 
the angular dependence of Fermi-surface cross-section area in external magnetic field are obtained 
for various typical crystal symmetries. The simple azimuth-angle dependence of the Yamaji angles 
is derived for the elliptic in-plane Fermi surface. These formulas correct some previous results and 
allow the simple and effective interpretation of the magnetic quantum oscillations data in cuprate 
high-temperature superconducting materials, in organic metals and other Q2D metals. The relation 
between the angular dependence of magnetoresistance and of Fermi-surface cross-section area is 
derived. The applicability region of all results obtained and of some previous widely used analytical 
results is investigated using the numerical calculations. 
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I. INTRODUCTION 

The layered quasi-two-dimensional (Q2D) compounds attract great attention for their novel physical properties and 
promising technical applications. High-temperature cuprate superconductors, [l| organic metals, heterostructures,[3j 
intercalated graphitesjj] are the examples of these compounds. The knowledge of quasiparticle dispersion in these 
compounds is very important for understanding their properties and electronic phase diagram. The traditional and 
powerful tools to determine the Fermi surface (FS) geometry and the electron dispersion in various metals are the 
magnetic quantum oscillations (MQO)@ and the angular dependence of magnetoresistance (ADMR)[5j. There is a 
huge amount of publications, devoted to the experimental determination of the FS geometry and electron dispersion 
in high-temperature cuprate superconductors 0,13: m MgB2,[9] in organic metals (see refs. (o|.[lo| for reviews) and in 
many other Q2D metals. The interpretation of the MQO data is, usually, based on the detailed comparison with the 
band-structure calculations, which is a complicated and often ambiguous procedure. The interpretation of ADMR is 
also based on fitting by the numerical calculations with a large number of fitting parameters, [l lMl 3j j The quick and 
effective extraction of the FS geometry and of electron dispersion from the experimental data on MQO and on ADMR 
requires reliable and simple theoretical formulas. 

The general form of the electron dispersion in Q2D compounds with monoclinic or higher crystal symmetry can 
be expressed as the Fourier series in cylindrical coordinates: 

e(k)= ^2 e f iU s ) cos {vkzC*) cos ([i<f> + 4>iiv) i (1) 

i/>0,/i— even 

where the integers v, fi > 0, the electron momentum k = (k x ^k y ,k z ) (we put h — 1), c* is the interlayer lattice 



constant, k — + k y is the absolute value of the in-plane momentum and </> = arctan (k y /k x ) is the in-plane angle 

(i.e., the azimuth angle in spherical coordinates). In triclinic crystals the only symmetry constraint on the electron 
dispersion is e (k) = e (— k), and the electron dispersion ([T]) may also contain the additional terms 

Ae(k) = ^ ^2 (fy sin (vk z c*) cos + MV ) . (2) 

f >0 fi— odd 

For simplicity, below we only consider the case of monoclinic or higher crystal symmetry, where the terms @ arc 
absent. 

Usually, it is sufficient to keep only the first few terms in the infinite series ([T]). For example, if the interlayer transfer 
integral of conducting electrons, t c (<j)) ks eoi, is much smaller than the in-plane band width eoo, the tight-binding 
approximation can be used, and one keeps only the terms with v = and v = 1: 

e(k) =e(M) -2M0)cos(fc 2 c*). (3) 

The FS, being given by the equation e (k) = Ep, is a warped cylinder in Q2D compounds. If magnetic field is applied 
along the z-axis of the Q2D metals with the electron dispersion ©, there are two extremal FS cross-section areas A ext 
encircled by the closed curves e {k x , k y )±2t c = Ep. Hence, the two close fundamental frequencies F± 2 = (c/2neh) A ext 
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appear in MQO, giving the beats of MQO.@ The temperature dependence of the MQO amplitude gives the cyclotron 
mass for the extremal orbit: m* ext = (l/2ir) [dA ext /dE]. The beat frequency AF = F\ — F 2 gives the interlayer 
transfer integral: At c = AF (eti/m* xt c). The difference of the two extremal cross-section areas AA ext = 1-neh AF j c 
and, hence, the beat frequency depend on the magnetic field direction. In the first order in the interlayer transfer 
integral and for the axially symmetric FS, this dependence is given bv|14j 

AA ex t oc Jo (c*kp tan 9) , (4) 

where Jo is the Bessel function, fcp is the in-plane Fermi momentum and 9 is the tilt angle of ma gne tic field B with 
respect to the z-axis (the polar angle of B). Eq. Q was first derived geometrically by Yamaji [14| to explain the 
oscillating angular behavior|15j of interlayer magnetoresistance in Q2D organic metals. As the difference between the 
two extremal cross-section areas is proportional to the interlayer transfer integral t c (9), Eq. @ suggests that the 
interlayer transfer integral has the similar angular dependence: 

t c (9) «t c (0) Jo (c* A* tan 0), (5) 

which gives a strong angular dependence of interlayer magnetoresistance p zz oc xr c (9). Eq. ([5]) was later confirmed 
by the quantum-mechanical calculation of the amplitude of interlayer electron tunnelling in tilted magnetic field. [l6| 
The angles 9 m , for which the Bessel function has zeros: 

Jo (c*kFtan9 m ) = 0, (6) 

are called the Yamaji angles and used to determine the in-plane Fermi momentum kp. At these angles both the 
interlayer magnetoresistance and the amplitude of MQO have maxima. Usually, the in-plane electron dispersion 
£{k x ,ky) is anisotropic, and Eqs. (J4j) , ([5j) acquire a (^-dependent correction, where if is the azimuthal angle of the 
magnetic field direction, tamp — B y /B x . There is a considerable practical need of the simple analytical formula 
for the ip- dependence of AMRO and MQO, which can be used to extract the in-plane electron dispersion from the 
experimental data. 

The widely used analytical result for the (^-dependence of the FS cross section, derived by Bergemann et al. [|[ and 
given by Eq. (I24[) below, takes the FS corrugation only in the first order, which is not enough to obtain correctly 
even the main (/j-dependent term in the angular dependence of the cross-section area. Another simple and widely 
usedp7l-[l9j analytical result for the (^-dependence of AMRO maxima (Yamaji angles), 

tan0„«^(n-l/4)/pg ax c*, (7) 

with p n ax being the maximum value of the Fermi momentum projection on the in-plane magnetic field direction, was 
derivedjl7j from the Shockley tube integral[Il| using the saddle point approximation. This approximation assumes 
that the z-component of the electron velocity oscillates rapidly when the electron moves along its closed classical orbit 
in the momentum space in magnetic field. This is valid only at high tilt angles 9 of magnetic field, when tan# ^> 
l/c*p^ ax (ip), and only in the very clean samples with uj c t cos 9 3> 1, where w c is the cyclotron frequency and r is the 
electron mean free time. For the first Yamaji angle this derivation is too approximate because c*p^ ax (ip) tan 9 ~ 1. 
Below we show that Eq. ([7]) is valid only for the elliptical FS in the limit WcTCOs^ 3> 1. With some small error 
it can also be applied to the FS, which is close to elliptical. However, Eq. ([7]) gives completely wrong result for 
the (^-dependence of Yamaji angles when the in-plane FS has tetragonal (as in cuprate high-Tc superconductors) or 
hexagonal (as in MgB2@ or intercalated graphites [H) symmetry. 

The aim of the present paper is to derive the new suitable analytical formulas for the ^-dependence of the FS 
cross-section area, Yamaji angles and magnetoresistance, which can be used to extract the FS parameters from the 
experimental data. The applicability region of some previous and widely used results will also be studied. 

In Sec. II we write down the relation between the dispersion (fTJ) and FS harmonic expansion. In Sec. Ill we find 
the main (^-dependent correction to the FS cross-section area for the anisotropic dispersion e (k x , k v ) in Eq. ([3]), when 
the in-plane anisotropy of the FS is weak. As will be shown, this result has wide applicability region and can also 
be applied to almost square-shaped in-plane FS as in the high-Tc cuprate superconductors. In Sec. IV we derive 
the exact expression for the Yamaji zeros for the elliptical FS shape. The deviations from this result for non-elliptic 
FS will also be studied. In Sec. V we derive the relation between the fc z -dependence of the cross-section area and 
magnetoresistance in the clean samples, where uj c tcos9 ^> 1. This relation shows, that the geometrical and resistivity 
Yamaji angles coincide in the limit uj c t cos 9 » 1. The discussion and summary of the results is given in Sec. VI. 
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II. FERMI SURFACE PARAMETRIZATION 



The dependence of the Fermi momentum kp ((f), k z ) on the polar angle <fi and the momentum component k z can be 
expanded in the Fourier series: Q 

k F (tf),k z ) = ^2 k v (</>) cos (vk z c*) (8) 
= ^ v cos (^^ C *) cos + • ( 9 ) 

Ai,i/>0 

The Fermi momentum satisfies the equation 

s[k F ((j),k z ),^k z ] = E F , (10) 

where Ep = e (kp) is the Fermi energy. The coefficients in Eq. (JlJ are related to the coefficients fc M „ in the FS 
parametrization © through the equation 

^ e M „ [k F (<j), k z )] cos (vk z c*) cos Qu^ + = Ep. (11) 

i/>0,/i— even 

This equation on k^ can be solved by the iteration procedure, assuming that the warping coefficients e^/soo are 
small and fall down rapidly with increasing fi and v. In the first order, each term k^ v in the series © comes only 
from the term in ([I]) with the same indices fi, v. 

k$ = -e^(k F )/e' 00 (k F ). (12) 

In the second order in e^„, the coefficients k^ v come from the interference of the infinite number of the terms e M v' 
and e^z/jy// in the dispersion ([T]), such that /i = jjl ± p!' and v = v 1 ±v" . 

For simplicity, we take the dispersion © and assume that the k z dependence of the energy is weak, i.e. the interlayer 
transfer integral t c <C Ep. The solution of equation 

e (k, 4>) — Ep + 2t c cos(k z c*), 

in the first order in the interlayer transfer integral t c = t c ((f>) , gives the FS shape in the cylindrical coordinates: 

k F (cf), k z ) = k ((f)) + ki ((j>) coa(k z c*), (13) 

where k (<fi) satisfies e [k ((f)) , <f)\ = Ep and 

h ((f)) = 2t c ((f,) I [de (k, (f)) /dk] \ k=koW - (14) 

The partial derivative (de/dk) \k=k {<f>) is the projection of the Fermi velocity on the line, connecting the point on the 
FS with the coordinate origin k = 0. It depends on the electron dispersion and on the azimuthal angle (f>. 
For the quite general form of the electron dispersion, 

e(fc,0) = fc Q 5 (0), (15) 

where g (<$>) is an arbitrary function and a is also arbitrary, the derivative 

[de (k, 0) /dk] \ k=k()W = Ep/k Q ((f)) . (16) 

The superelliptic dispersion 

e(k x ,k v ) = (k x /k 1 ) a + (k y /k 2 ) a , (17) 



which includes both linear and quadratic dispersions, is only a particular case of the dispersion (j 15[) . With the relation 
(EUD , Eq. Jin]) simplifies to 



k F (4>,k z ) = k ((f)) 



l + ^Mcosfec*) 

Hip 



(18) 



However, the relation (|16p may violate in some compounds, and the application of the simplified formula (|18[) instead 
of Eqs. (IT3"l) . (fH| requires additional proof. 
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III. CROSS-SECTION AREA 



If the magnetic field is applied at polar and azimuthal angles 9 and tp, the Fermi surface cross-sectional area 
A = A {k Z Q, 9, tp), cutting the z-axis at k z = k z g and perpendicular to the field, is given by the integral 

A(k z0 ,9,^) = #' kp (ip + <p',k z )/ (2 cos 0), (19) 



where <p' is the angle in the x-y plane between the direction of magnetic field and a point on the FS, and k z at this 
FS point satisfies the equation 

k z = k Z Q — kp (ip + (/)', k z ) tan 9 cos </>'. (20) 

Eqs. (fl~9| and (|20| allow to find the cross-section area A (k z o,9, ip, kpc*) numerically for any given FS, determined 
by the function kF (</>, k z ) or, equivalently, by the coefficients k^ in the expansion (|9|). In practice, one usually solves 
the inverse problem of the extraction of FS parameters from the experimental data on MQO or AMRO. Then, the 
direct procedure of fitting the experimental data by the parameters in the expansion (J9j) is rather ambiguous because 
of too large number of fitting parameters. Usually, the coefficients k^ v fall down rapidly with increasing /i and v. 
Therefore, it is useful to fit only the first few terms in the similar harmonic expansion of the cross-section area 

A (k z0 , 9, ip) = ^2 A iiv (#) co s [W + M cos (vc*k zQ ) , (21) 

keeping only the first few terms in the expansion ©. The first coefficients A^ (9) can be found analytically in the 
main order in k^ v . The analytical formula for the coefficients A^ (9) of the cross-section area is especially useful 
because of their rather complicated dependence on 9. In Sec. V it will be shown that the coefficient A^i (9) is directly 
related to the angular dependence of magnetoresistance at lu c t ^ 1 and t c /Ep <C 1. 

In the zeroth order in coefficients k^ v in the expansion ([9]), i.e. for cylindrical FS neglecting any warping and 
in-plane asymmetry, one obtains the trivial result A^ — nkp/ cos 9, where fci? = fcoo- In the first order in these 
coefficients k fiu% one can neglect the dependence kp (</>, k z ) in Eq. (I2U1) and substitute k z — fc z0 — kp tan 9 cos </>' to Eq. 
©. Then, substituting Eq. © to Eq. (IT9l . one obtains the first order correction to A^>: 

J 2 cos 9 

27r fcoo<#' r i i a'\ , a l 

CQsd 2^ k^cos[n{tp + (j)) + cf) fl } 

H,v>Q 

x cos \v (k z0 — kp tan 9 cos (j) 1 ) c*] , 

where the sum v>n ' does not include the term \i = v = 0. Since [i is even, one can replace in the integrand (here 
and later we introduce the notation k = kpc* tan#) 

cos [v (k z o — kp tan 9 cos cj>') c*] — > cos [^fc z oc*] cos [vk cos <fi'] . (22) 

One can also replace in the integrand 

cos [fi (ip + (j)') + cos [fnp + M ] cos \p4>'\ , (23) 

because all odd terms vanish after the integration over (j)' . Then, after the integration over </>', the correction A^ in 
the first-order in k^ v writes down as 

A(1) = ^ ,(_1)2M C ° S [w + ^ C ° S {l/kz0C * ) ^ (24) 

in agreement with Eq. (2) of Ref. [8[. Since J M (0) = for /i ^ 0, all terms ~ fc^o vanish in (|24|) . This is natural, 
because in the zeroth order in t c the cross-section area 

A^(k z0 Av)= I'" d<t>' kl <? + P (25) 
J 2 cos 9 



5 



is independent of p. Hence, to extract any information about the (^-dependence of the FS, one needs to consider the 
first order in t c /Ep, i.e. to find (8). Thus, the ^-dependence of the cross-section area starts from the term kui, 
which is of the same order as the second order term k^koi/kp [see Eq. ([TBjl], Since Eq. (f24|) . or Eq. (2) in Ref. 
is derived only in the first order in k^ v , it does not give the correct (^-dependence of the cross-section area even in 
the lowest (^-dependent order. This is illustrated below in Figs. II 121 The extraction of the higher harmonics using 
Eq. (|24|) is even more incorrect. 

Let us calculate more accurately the lowest-order (^-dependent term in the cross-section area, which is given by 
the coefficient A^i (8) in the Fourier expansion. To calculate this coefficient in the main order in FS warping, it is 
sufficient to use the FS shape in the first order in t c , given by Eq. (|f 31) . Then, in the same order, Eq. (j!9p rewrites 



A(k z0 ,9,(p) 



2cos6> 



k (ip + <p') 



(26) 



and substituting Eqs. (|2T)|) and (|22l) . we obtain the following expression for correction to A (k z o, 9, <f)'- 



cos( 

Here we have also changed the integration variable: 0' — > 0' + tp. The Yamaji formula (J4j) is easily obtained from 
(|2T|) after taking ko (0) — kp — const, fci (0) = C\ (2t c /Ep) kp = const, where the dispersion-dependent constant 

d = (Ep/kp) I (deoo/dk) \ k=kp ~ 1, (28) 

and the integration over 0', resulting to 

The lowest-order (^-dependence of the cross-section area is determined by the Fourier coefficient A m \ (8), given by 



2tt 



COS 



(mip + ip ml ) dip A^l (k z0 , 8, ip) 



A m i (8) = / " v " ' " ' ■ (30) 

J 7rcos6> cos [c*k z0 \ 

Performing the integration over ip, we obtain 

Am (8) = K ' / dc/>'k (00 h (00 cos (m0' + <p ml ) J m [c*fc (00 tan 8} . (31) 
cos 8 J 

To go further, we need to specify the functions fco (0) an d fci (0). We distinguish two symmetries of electron dispersion, 
namely, with straight and (^-dependent (corrugated in the main order) interlayer transfer integral. 



A. Straight interlayer hopping 

When the in-plane FS anisotropy is weak, one can keep only the first 0-dependent term in the Fourier expansion of 
the functions fco (0) and k± (0). If the crystal symmetry allows the 0-independent (straight) interlayer coupling, these 
functions expand as 

fc o (0) ~ (1 + /3cosm0)fc F , (32) 
It C 

fci (0) ~ — f — (1 + Pi COS 7710) kp, 
Ep 

where m is an even integer number, , <C 1 and the constant C% ~ 1 is given by Eq. (f28|) . Now we expand 
J m [c*fco (00 tariff] in the small parameter f3 up to the first order (for the first Yamaji angle k = c*kp tand » 2.4 ~ 1), 
and the integral over 0' in Eq. (|3"Tj) simplifies to 

dcf)' (j3 + Pi) cos (m0O cos (m0O J m («0 

n 

d<j>' cos (m0O J' m (k) k/3 cos (m0O 
vr/3 { J m (k) (1 + [ii/P) + J' m (k) k} , 
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where the derivative 



ji t \ dJ m (k) m . . . , 

J m { K ) = = —Jm («) - 4+1 («) 



Hence, in the first order in /3 we obtain 

Anl (^) = (-I)'' 



l/2 4:7rk 2 F Ci/3t c 



Jm (k) 1 + ^ + 



/3 



-Ef cos 

Combining the results (j2"Tj) . (|29l) and (f3"3")l . we obtain the cross-section area 



m+1 



A(fc zO ,0,¥>) 



COS -Bp COS ( 



cos [c*/c z0 ] 



{ Jo (k) + f3 (-l)" l/2 [(1 + + m) J m (k) - «J m+1 («)] cos (my)} 



(33) 



(34) 



The ratio Pi/ ft, entering this formula, can also be expressed via the FS parametrization, given by Eq. j3 = 

fc M o/fcoOj Pi — fejui/fcoii Pi/P = k^ikoo/koik^Q. The constant Ci is equivalent to the renormalization of t c and does 
not influence the Yamaji angles. However, it changes the beat frequency of the magnetic quantum oscillations. It 
also changes the amplitude of the (/^-dependent term in the cross-section area. For the dispersion of the form (|15p . 
i.e. if the relation P^|) satisfies, the ratio /?i//3 = 1 and the constant C\ = 1. For arbitrary dispersion, Pi/j3 ~ 1 and 
d ~ 1. 
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FIG. 1: (Color online) The in-plane FS shape (a) and the first Yamaji angle as function of the azimuth angle <p (b) for the 
elongated FS with monoclinic symmetry and the Fermi momentum given by ko (0) = 2~^ 4 =o cos (4?^) with f3 = kio /koo = 
0.5. The higher harmonics are added to smooth the FS (see Fig. a). In Fig. b the blue solid curve is the numerical result for 
the first Yamaji angle, obtained using Eq. ()27[) ; the red dotted curve is the result of Eq. (I45|l ; the green dashed curve is the 
result from Eq. (|34[) and magenta dash-dotted curve is the result form Eq. (|24[) derived by Bergemann et al.[8fl One can see 
that Eq. (|45p gives the best result for Yamaji angles when the FS is strongly elongated. The harmonic expansion still gives 
a reasonable result, though the coefficient /3 = 0.5 ~ 1 is not small. The result of Eq. (|24[) gives much weaker (^-dependence, 
than the correct one. 



The difference between two analytical results, given by Eqs. (f3"4"|) and (|24l) . is very strong: the factor J m (k) in 
Eq. (IM1) is replaced by the completely different factor [(1 + (3i/(3 + m) J m (k) — nJ m +i («)] in Eq. ([M]). First, the 
^-dependence of the cross-section area, predicted by Eq. fl34[) . is stronger approximately by a factor 2 + m than 
that of Ref. Q. Second, it may have different 0-dependence due to the J m +i (k) term, especially for high tilt angles 
k > 1. To illustrate the above statement, we plot the results of Eqs. (JTSJ), (|53)l. and (jtSJ) in Figs. Q] and ED for the 
dispersions with monoclinic m — 2 and tetragonal m = 4 symmetries. 



e Yara [Rad] 




FIG. 2: (Color online) The in-plane FS shape (a) and the first Yamaji angle as function of the polar angle (j> (b) for the FS 
with tetragonal symmetry and the Fermi momentum given by Eq. (|32[l with ft = —0.06. Fig. a shows, that even for very small 
parameter j3 the FS is strongly different from the cylinder. It is almost quadratic. Hence, keeping only the first term in the 
harmonic expansion of the ^-dependence and the assumption that this term is small, is applicable for most compounds with 
tetragonal symmetry. In Fig. b the blue solid curve is the numerical result for the first Yamaji angle, obtained using Eq. (|27|l . 
The red dashed curve is the result of Eq. (|45[) ; it gives too strong and opposite ^-dependence, which is completely incorrect. 
The green dotted curve is the result from Eq. I)34[): it gives very good agreement with the numerical result. The magenta 
dash-dotted curve illustrates Eq. ((24}, used by Bergemann et al.Q; this result gives too weak (^-dependence in agreement with 
the discussion in the end of Sec. Ill A. 



B. Strongly (^-dependent interlayer hopping 

In some compounds, e.g. in the high-temperature superconductors Sr2Ru04 and Tl2Ba2Cu06+5,[l, [Tl|-[l3| the 
body-centered tetragonal symmetry of the crystal leads to the (^-dependent lowest-order interlayer transfer integral, 
t c (4>) = t c o sin (2(f)), in the all or some parts of the FS. Then, instead of Eq. (|32|) . we have 

k (</>) « (1 + /3 cos 2TO0) k F , (35) 
2t 

ki (cf)) s» — ^-kpCi sin (m4>) (1 + (3\ cos2m</>) . 
Ef 

Substituting this into Eq. (|3"Tj) we obtain the main ^-dependent term, determined by the coefficient 

Aml {9) " — — W Jm (K) (36) 
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in agreement with the first-order result, given by Eq. (|24[) . This term does not depend on the in-plane FS anisotropy 
(3. To extract this anisotropy in the first order in the (^-dependent interlayer transfer integral, one needs to consider 
3m harmonic in the cross-section area. For this we replace in Eq. (|3Tj) m 3m, substitute Eq. (|35j) and perform the 
calculation, similar to that in the derivation of Eq. (1331) . Then, in the lowest order in /3,/3\ we obtain 



A 



3m 1 



(0) = 



(_1) 3 "V 2 2 Tit c C 1 
cos 9 E F 



1 + — + 3m J 3m (re) - k J 3m+ i (re) 



(37) 



This result differs from Eq. (f33|) by the replacement m — >• 3m (note, that m = 2 for Sr2Ru04 and T^E^CuOg+a), 
and the prefactor before the square brackets is two times smaller. The difference between the first-order result of 
Eq. and the new formula (|3"7|) for the 3m harmonic is even stronger than in the case of Eq. (|3"3")l . The total 

(^-dependence of the cross-section area in the case of (^-dependent interlayer coupling, given by Eq. (|35p , writes down 
as 



A(k z0 ,d,ip) p 
x<j J m ( K )sin(m^) + ^(-l) 3m/2 



4nk F t c Ci 



cos -Bp cos 

/3 



1 



3m 



cos [c*/c z o] x 
1/ 

^3m (k) - KJ 3m+ i (re) 



(38) 



sin (3m<£>) 



This formula can be applied to analyze the experimental data in high-temperature superconductors Sr2Ru04, 
T^E^CuOg+a, where m = 2 in Eq. (|38[) . and to some other layered compounds with the appropriate symme- 
try. Note, that Eqs. (j34)) and (|38)l were derived under condition (3k <C 1, which is fullfilled in the compounds of 
tetragonal or hexagonal symmetry at not very high tilt angle of magnetic field. At very high tilt angle, tan# > Ep/t c , 
above derivations are not valid also because of the multiple intersections of the FS by the cross-section plane. 



C. Analysis of magnetic quantum oscillations 



The well-resolved magnetic quantum oscillations in quasi-2D metals give two close frequencies -F max and F m i n , 
corresponding to the maximum and minimum of the FS cross-section area. To extract the (^-dependence of the FS, 
as follows from Eqs. and (f3"5|). one needs to measure the (^-dependence of the difference AF — F max — F m j n ~ 
Ankptc/ Ep between these two close frequencies (the beat frequency), which is harder because requires the resolution 
of MQO in the wider interval of magnetic field. The observation of the beat frequency itself is important, because it 
means the existence of the 3D Fermi surface, i.e. of the coherent interlayer electron transport. Eqs. (|3~4"f and (|3"5)) 
can be used to determine the optimal orientation of magnetic field for the observation of the beat frequency. For 
straight interlayer electron coupling as in Eq. (1321) . the beat frequency has maximum value when magnetic field is 
perpendicular to the layers, i.e. at polar angle 9 = 0. However, for the ^-dependent interlayer coupling at 9 = 
the beat frequency is zero, as follows from Eq. ([55]) . Hence, in this case to observe the beat frequency one needs 
to incline the magnetic field. The angular dependence of the beat frequency is given by the function in the curly 
brackets in Eq. (|38p . The first term J m (re) sin (mip) in the curly brackets is much larger than the second. Its 
maximum gives the optimal orientation (9 op t, f pt) of magnetic field for the observation of MQO beat frequency. For 
m = 2, as in Sr2Ru04, T^E^CuOg+s and some other high-Tc compounds, the factor J m (re) sin (imp) has maximum 
at tpopt — (2ti + l)7r/4 and 9 op t ~ arctan (3.0/kpc*), where c* is the interlayer lattice constant. Note, that the spin 
factor of MQO also depends on the angle 

If the beat frequency of MQO cannot be resolved (in dirty materials or at high temperature), the minima of the 
beat frequency, i.e. the Yamaji angles, can be detected from the increase of the amplitude of MQO. This increase 
of MQO amplitude happens because at the Yamaji angles the MQO from both extremal electron orbits have the 
same phase. 0, Q The Yamaji angles can be much easier distinguished from the angular dependence of background 
magnetoresistance (AMRO) (see Sec. IV). To determine the (^-dependence of the Yamaji angles one can again use 
Eqs. (JM]) and (j55]>. 

If the (^-dependence of MQO beat frequency is clearly resolved, one can obtain the information about the in-plane 
FS. Eqs. ((34]) and (|38[) again can be used to determine the optimal magnetic field orientation for the observation of 
this ip- dependence. In the case of straight interlayer electron hopping, this (^-dependence cx cos mip has the maximum 
amplitude when the factor / (re) = (1 + /3i//3 + m) J m (re) — ftJ m +i (re) in Eq. (|34p has maximum. For typical value 
/3i//3 = 1 this factor as function of re = kpc* tan# for m = 4 and m = 6 is plotted in Fig. [3] The function / (re) has 
first maximum at re ~ m for m = 4 and m — 6 (see Fig. [3]). It is reasonable to use only the first maximum, because 
at high tilt angle of magnetic field the cyclotron mass is large and the amplitude of MQO is too small. 

In the case of straight hopping, already the lowest-order harmonic in the (^-dependence of the MQO frequency 
gives the relative amplitude /3 of the same harmonic in the ^-dependence of the in-plane Fermi momentum [see 
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FIG. 3: (Color online) The amplitude of the (^-dependence of the extremal cross-section area as function of the polar angle 9 
for tetragonal (blue solid line) and hexagonal (green dashed line) symmetries in the case of straight interlayer electron hopping. 



Eq. (|34l) ]. In the case of ^-dependent interlayer electron hopping, given by Eq. (|35p . in the main order, the in- 
dependence of MQO frequency comes from the ^-dependence of the interlayer transfer integral and does not give 
information about the in-plane FS. To determine the shape of the in-plane FS, one needs to study higher harmonics 
oc cos (3m<p) in the MQO frequency. The amplitude of the cos (3m<p) term in MQO frequency is given by the function 
/i (k) = (1 + f3i/(3 + 3m) J^ m (k) — KJ3 m +i (k) in Eq. (|38|) and has first maximum at k w 3m. This determines the 
optimal polar angle 9 optl at which this ^-dependence is most easily observed. According to Eq. (|38[) . this dependence 
gives the amplitude of the 2m harmonic modulation of the in-plane FS. 



IV. ELLIPTIC FERMI SURFACE 



Now we derive the analytical formula for the Yamaji angles for the elliptic in-plane dispersion 

e(k x ,k y ) = k 2 x /2m x + k 2 /2m y =e(k) [1 + /3cos20] , (39) 

where e (k) = k 2 /2m, m = 2m x m y / (m x + m y ) and /3 = (m x — m y ) J (m x + m y ). The shape of the FS for this 
dispersion is, of cause, also elliptical. The ellipse can be obtained from the circle by the dilation A^ along one in-plane 
direction (along the x-axis): x — > Ax. Consider the cross-section area of the FS by the plane, cutting the fc z -axis at 
the point k z $, and perpendicular to the magnetic field direction B = Bn, where the unit vector 

n = (n x ,n y ,n z ) = (sin 9 cos (p, sin 9 sin ip, cos 9) . (40) 

For the circular in-plane FS this cross-section area is independent of the angle ip. In the first order in t c / Ep, it is 
also independent of k z o at special directions ny am , corresponding to the Yamaji angles 9 — 9Yarm given by Eq. ([6]). 
After the dilation K x , the direction of magnetic field, which is perpendicular to the cross-section plane, also changes: 

n ^ Al („) = ^teL. (41) 



y/(n x /\) 2 +nl+: 



However, the cross-section area perpendicular to ni = A x (ny am ) remains independent of fc z o, if it was independent 
before the dilation. Hence, the direction ni = A^ (ny am ) corresponds to the new Yamaji angle 0y am ((p). The polar 
and azimuthal angles are related to the components of the vector ni = (ni Xl n\ v , n\ z ) as 



t&n9 1 = - , tan</>i = — (42) 

nu ni x 

Combining above equations we obtain the relation between the old and new Yamaji angle 9y- am — A x (#y a m) 



tan^ am = V^/ A2 + j = /cosV | ~ 
ta,n 9Yam n z tan 9y a m V A 2 
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The angle <p> here is the angle before the dilation A^. It is related to the angle <pi after the dilation as 

tan ipi = A tan ip. 

Then, after simple trigonometric algebra, we obtain 

tan^ nm _ cos^ r | — ~ _ _ 1 

tan VYam A V ^ 2 + tan2 V ^ 2 cos2 Vl + sm2 

For the elliptic dispersion (|39p the maximum value of the Fermi momentum projection on the in-plane magnetic field 
direction is given by 

PT X = Vbi co ^) 2 + (P2sin^) 2 , (44) 

where p\ = 2m x EF and p| = 2m y £F- The r.h.s. of Eq. (|43|) coincides with P2/p^g ax - Hence, the generalization of the 
Yamaji zeros to the elliptic dispersion (f3l)| writes down as 

J o [c*pg ax (^)tan0 n ]=O. (45) 

Approximately, Eq. (|45|) coincides with Eq. ([7]), derived for the interlayer conductivity^ [13] from the Shockley 
tube integral[2l|. The saddle point approximation, used in Ref. [13] to derive Eq. ([7]), assumes that the z-component 
of electron velocity oscillates rapidly when the electron moves along its closed orbit in the momentum space. This 
is valid only at high tilt angles 9 of magnetic field, when tan 9 3> l/c*p'^ ax (<p), and only in very clean samples with 
WcTCOs^ S> 1. The reason, why Eq. ([7J describes well some experimental data,[f| comes from its coincidence with 
the exact geometrical expression (|4"5j) for the Yamaji angles for the elliptic Fermi surface, which according to Eq. (|55|) 
gives the maxima of magnetoresistance. 



V. MAGNETORESISTANCE 



To calculate magnetoresistance as function of the direction (9, ip) of magnetic field one can use the quasi-classical 
Boltzmann transport equation for electrons moving along the closed orbits in magnetic field. This approach gives the 
Shockley-Chambers formula pi), which at zero temperature expresses conductivity tensor o a p via the integral over 
the Fermi surface: 

, n * e 2 f „ m% cos# /u>h 

/■2tt p2-7T 
JQ JO 

Here the momentum space is parametrized by the momentum component kn — fc z o cos 9 along the magnetic field, by 
energy E and by the angle < ip < 2ir of the rotation in the cross-section plane. The effective cyclotron mass of the 
orbit is given by 

* = J_^i - 1 OA 
m " ~ 2vr 8E ~ 2ne' (k F ) dk F ' [ } 

where A — A (kjj , Ep) is the area of the FS cross-section perpendicular to the magnetic field B at the momentum 
kii || B. The cyclotron frequency of the orbit u>h = eB /m* H c, and v a (ip, ku) is the component of the electron velocity 
on the FS. Generally, the mean scattering time r in the integrand (|4l)|) may also depend on the position on the Fermi 
surface. However, we neglect this dependence because in the simplest theory of spin- independent short-range impurity 
scattering r depends only on the density of states at the Fermi level. 

For dispersion ([3]) the electron velocity component along the z-axis is a function of k z only: 

v z {k z ) = (2c*t z /H)sm(c*k z ). (48) 

The k z coordinate of the FS point K satisfies the equation (|20|) . where cp + <f)' is the azimuthal angle of the projection 
of the FS point K on the x-y plane. Approximately, this equation can be solved by the iteration procedure. In the 
zeroth order 

fc(°) = k z0 - k F cos (</>') tan 9, (49) 
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and in the next orders 

/c< 4+1 ) = k zQ - k F (V + <p, fcW) cos ((/)') tan d. (50) 

The angle ip entering the Shockley-Chambers formula P5|) corresponds to the increment of the cross-section area 
at a given increment of energy: 

1 dk dk dk± 
d v = — = ; o F - 

m H v± m H oE 

Generally, ip is different from the angles (p' and (p of the rotation in the cross-section and in the x-y planes. The 
cross-section area multiplied by cos 9 is equal to the area of the projection in the x-y plane, and ip is related to the 
angle <p of the rotation in the x-y plane as 

dip _ k F ((f), k z ) dk F Q, E) 

d(f> ~ m* H cos9 BE ' [ ° ' 

For cylindrical FS one has k F (4>) — k F , and ip coincides with the angle (p. 

Now we show that for ujht 3> 1 the minima of a zz (9,ip), given by Eq. (|46|) . coincide with the minima of the 
mean-square value of the derivative (dA/dk z o), i.e. with the geometrical Yamaji angles. At loht 3> 1 the exponent 
e.-VI>*iiT ~ i : an d Eq. flU gives 

<Taa(9,<p) = A^I dk z0 - m H C ° Sd /"" ( f\ a kzQ ) d S . (52 ) 



An^h 2 J 1 — exp (—2tt/u)ht) 
Using Eq. (fBTj) we transform the integral 

2-7T 

dipv z (ip, k z0 ) 



o 

m* H cos 9 dE dk z 

2T (</>, k z ) dk F (</>, k z ) 

m* H cos 9 dk z 

The derivative 

dk F (0, k z ) _ dk F [cp, k z (k z0 ,(p)} 
dk z dk z0 ■ (dk z /dk z0 ) ' 

From (PH)) in the first order in t z we obtain dk z /dk z o = 1. Hence, from (I5"3")l we get 

r2 " dk F (cP,k z ) _dA(k z0 ,9, V o) 



(53) 



(54) 



I = 



m* H cos 9 2dk z o dk z o m* H 



where the cross-section area A is given by Eq. ([T 
Now from (|52|) we get 



e 2 rcos6» f dk z0 f dA(k z0 ,9,ip)\ 2 

'» {9 > = "w J ^ { dkz0 ) • (55) 

Similar relation without rigorous proof was also proposed in Ref. [l7|- Eq. ([55)1 means, that the angular dependence 
of the intcrlaycr conductivity a zz and of the mean-squared derivative of the FS cross-section area dA/dk z o coincide in 
the limit ujht ^> 1. In particular, the geometrical Yamaji angles coincide with the minima of interlayer conductivity 
at ujht 1. 



12 



VI. SUMMARY AND DISCUSSION 

Above we have obtained the following main results: (i) the analytical formulas (|34|) , (|38|) for the main (^-dependent 
term of the cross-section area A (fc 2 o, 9, <p), when the FS corrugation is weak (Sec. Ill); (ii) the exact analytical formula 
(|7|) for the Yamaji angles in the case of elliptic in-plane Fermi surface (Sec. IV); (hi) the derivation of Eq. ([55]) . which 
states that in the limit lort 3> 1 the angular oscillations of magnetoresistance coincide with the angular oscillations 
of the ^-dependent term in the cross-section area. Eq. ([55]) brings additional importance to the results in Sees. Ill 
and IV for the FS cross-section area. In particular, Eq. ([55]) means, that the geometrical Yamaji angles dy a m (p) 
coincide with the maxima of magnetoresistance at ujht ^> 1. 

If the interlayer electron hopping is </>- dependent, as in Eq. ([35]) , Eq. ([55]) also suggests the very strong in- 
dependence of the interlayer magnetoresistance, given by Eq. ([38]). Note, that this (^-dependence of magnetoresistance 
for Sr2Ru04 and Tl2Ba2Cu06+<5 is in contrast to the so-called "third angular effect", developed in Ref. [25| for the 
in-plane magnetic field direction. For the Fermi surface in Fig. the third angular effect [25] predicts maxima of 
conductivity at tp = 0, coinciding with the positions of the FS inflection points, while Eq. (|38|) predicts these maxima 
at (p = 7r/4. 

The 0—ip dependence of the cross-section area is also an important result, because the magnetic quantum oscillations 
give the extremal cross-section area as function of magnetic field orientation. In the end of Sec III we summarized, 
how the formulas (f34|) and (l38l) for the cross-section area can be applied to analyze MQO. In particular, these formulas 
give (i) the optimal direction of magnetic field for the observation of the beats of MQO and of the (^-dependence of 
the beat frequency; (ii) the Yamaji angles, where the magnetoresistance and the amplitude of MQO have maxima; 
(iii) the relation between the FS shape and the ^-dependence of the beat frequency. 

Now we discuss in more details the applicability region of all above results and compare them with the previous 
theoretical results. Let us first study the analytical formula for the cross-section area, obtained from the Fourier 
expansion in the limit of weak FS warping. This formula strongly differs from the previous result [8| [see Eq. (pM)) ] . The 
derivation of the result of Bcrgemann et al. 8] does not include the second-order terms in the FS corrugation k^/koo, 
which is necessary, because for the straight interlayer electron hopping and in the first order in the FS corrugation, 
the (^-dependence of the FS cross section disappears. The differences between the result of Ref. @ and the new 
formula are illustrated in Figs. 11121 Roughly, the formula of Bergemann et al.Q gives too small ^-dependence of 
the cross-section area, which is smaller than the correct result, approximately, by a factor of 2 + m, where m is the 
harmonic index of this ^-dependence. 

Eqs. (|34[) and (|38[) . derived using the small harmonic expansion, turns out to be a very good approximation for any 
typical FS with the tetragonal or hexagonal symmetry (see Fig. [3]) . This is because even a small value of j3 leads to 
the strong change of the in-plane FS shape: for the almost quadratic in-plane FS as in Fig. [2^, \(3\ — 0.06 <C 1. Eq. 
(|34[) may considerably violate only in the case of monoclinic or triclinic symmetry, when the in-plane FS is strongly 
elongated (see Fig. QJi). Fortunately, just for this case Eq. (pp)) gives the reliable result (see Fig. [T}d). 

The 8-ip dependence of the cross-section area, given by Eqs. (|34l) and ([3"5]) , allows to extract not only the leading in- 
plane (^-dependence of the FS, but also to get some information about the electron Fermi velocity from the coefficient 
Pi. However, the extraction of j3\ is much more difficult than the extraction of /3, because, according to Eqs. ([M]) 
and ([3"8"| . the dependence of the cross-section area A\ (6,(p) on fi\ is much weaker than on fi. In Figs. [TbUb, for 
definiteness, we take the dispersion of the form (fT5"|) and t c (</)) to be independent on <j>, which gives Eq. (|18[) and 
0i/P = 1, Ci = 1 in Eqs. (JSJ) and (38]). 

The (^-dependence of the Yamaji angles turns out to be weak in the cases of tetragonal or hexagonal symmetry. 
This dependence is much weaker than the prediction of Eq. (US")) , used in Refs. [I3-[2C| for the elongated FS, and is 
much stronger than the prediction of Eq. ([M]) . used in Refs. [8[. For the first Yamaji angle in the case of tetragonal 
symmetry, Eq. ([45]) even gives the opposite sign of the (^-dependence of the first Yamaji angle, as can be seeing from 
Fig. For the superelliptic FS, given by Eq. ([T7]) and discussed in Ref. [Hj], there is a strong tetragonal modulation 
of the elliptic FS, and Eq. (j^5]) also fails to give a reliable result. 

Eqs. (|34]) . ([38]) . (|45]) . (|55|) are valid only in the first order in the small parameter t^jep. This is, usually, a good 
approximation for the layered high-Tc superconductors, organic metals, and many other compounds. However, some 
fine details of the angular dependence of magnetoresistance may be sensitive to the next interlayer hopping term, 
especially in the case, when the main interlayer hopping is strongly ^-dependent. 

Eqs. (|34 p . ([38]) . (|45]) determine the geometrical conditions for the FS cross section to be almost independent on the 
interplane momentum k z , which, according to Eq. (|55l) . gives the minima of conductivity at uj c t — > oo. To check 
how strong the conductivity, given by the Shockley-Chambers formula ([4S]) . differs from the geometrical formula ([55]) 
at finite uj c t, we perform the numerical calculation of the ^-dependence of conductivity, a zz (#, tp, w c r), given by Eq. 
([46]) . at four different values of lo c t — 1, 2, 4, 8 and for two values of the azimuthal angle, ip = and tp = 7r/4. Then, 
we compare it with the dependence a zz (0,ip) given by Eq. (|55[) . The results for the polar-angle dependence of the 
normalized interlayer conductivity are given in Figs. 0] In this calculation we take the in-plane dispersion, proposed 
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for the layered cuprate high-Tc superconductors and given bv [lll - ll3j 

£ (k x , k y ) = 2ti [cos (k x a) + cos (k y a)] + 4i2 cos (^o) cos (k y a) — 2t% [cos (2k x a) + cos (2k y a)] — Ep, (56) 

where a = 3.95A is the lattice constant, t\ = 0.38eV, = 0.32fi, £3 = 0.5^2- We take the doping-dependent Fermi 
energy Ep = 0.02123ey. The FS for this dispersion (|56|) has tetragonal symmetry and is very similar to that in 
Fig. [2k with slightly different values kpa « 2.14 at ip = and kpa « 2.43 a,t (p = 7r/4. Therefore, we don't plot it 
again. The interlayer hopping term of electron dispersion is taken to be straight (^-independent) and given by Eq. 
([3"2"|) . From Fig. H]we see, that the geometrical formula ([53]) gives rather accurate results for the Yamaji angles of 
magnetoresistance at uj c t > 2. For large values of uj c t, the results of Eqs. (f55"1) and (|46|) coincides. At w c r < 1 the 
difference between Eqs. (|5"5")) and (|46l) for the first Yamaji angle reaches 5-10%. The (^-dependence of the Yamaji 
angles given by Eqs. (|55|) and (|46|) agrees very well; only the amplitude of AMRO reduces with decreasing uj c t. 




FIG. 4: The normalized conductivity u zz , calculated from Eq. (|46p for the dispersion (|56[1 as function of the polar angle 8 for 
two different azimuthal angles <j> at several values of uj c t = 1 (black dashed line), uj c t = 2 (red dash-dotted line), lo c t = 4 (blue 
short-dashed line) and ui c t — 8 (green dotted line). The solid magenta line gives the result of Eq. (I55[) . 

The background magnetoresistance and, in particular, the saturation values of <j zz at 9 — > 90°, depend strongly on 
the azimuthal angle if (see Fig. Therefore, it is reasonable to use this (^-dependence of the conductivity saturation 
value at 6 — ir/2 to determine the in-plane Fermi surface from the experimental data on the angular dependence of 
magnetoresistance. The theoretical prediction for this dependence can be obtained from the numerical calculation 
using the Shockle y-C hambers formula (|4"6")l , as is done in Fig. H) The origin of this (^-dependence is qualitatively 
explained in Ref. [25]. 
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